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In a recent paper [J, we investigated the effects of the massive charged gauge bosons (bileptons) 
predicted by the minimal 331 model on the off-shell vertex WWV* (V = 7, Z) using a SUl(2) x 
Uy{X) covariant gauge-fixing term for the bileptons. We proceed along the same lines and calculate 
the effects of the gauge bosons predicted by the 331 model with right-handed neutrinos. It is found 
that the bilepton effects on the WWV* vertex are of the same order of magnitude than those 
arising from the SM and several of its extensions, provided that the bilepton mass is of the order 

C ' of a few hundred of GeVs. For heavier bileptons, their effects on the WWV* vertex are negligible. 

The behavior of the form factors at high energies is also discussed as it is a reflect of the gauge 
invariance and gauge independence of the WWV* Green function obtained via our quantization 

(N ■ method. 

h : 

Ctf PACS numbers: 12.60.Cn,14.70.Hp,13.40.Gp 

S: 

r— I ' I. INTRODUCTION 

cm : 

, Radiative corrections to the WWV (V — 7, Z) vertex have long received considerable attention. Apart from its 
sensitivity to new physics effects, this vertex has theoretical interest by its own as it may serve as a probe of the gauge 
I \ sector of the standard model (SM). In this context, the one-loop contributions to the WW*y vertex, which defines the 
static electromagnetic properties of the W boson have been calculated in the SM 0, Q and several of is extensions 
such as two-Higgs doublet models (THDMs) 0], supersymmetric models @, 331 models @, 0], etc. Similar attention 
has been paid to the study of the WWZ vertex. Despite the attention has focused mainly on the static properties 
of the W boson, a more comprehensive study of the W boson properties requires the analysis of the off-shell WWV 
vertex, particularly when the neutral V boson is off-shell and the W bosons are on-shell as the resulting vertex could 
be tested with high precision at the CERN large hadron collider (LHC) or the planned future linear colliders. It is well 
known that although the on-shell WWV vertex renders a gauge independent amplitude by itself, gauge independence 
is lost once any external particle goes off-shell, thereby requiring the use of a nonconventional calculation scheme. 
Q_i The difficulties to obtain a gauge independent Green function for the WWV* vertex have long been known. In 
particular, a careless calculation of the radiative corrections to the WWV* vertex via a conventional gauge-fixing 
procedure irremediably leads to an ill-behaved gauge dependent Green function. For instance, one-loop corrections to 
the WWV* vertex were first calculated in the SM by the authors of Ref. [3J] using the Feynman-'t Hooft gauge via a 
conventional gauge-fixing scheme. The result was an infrared divergent Green function that violates unitarity. Aware 
of the fact that an off-shell Green function is generally gauge dependent, authors of Ref. [8[ invoked a diagrammatic 
■ method known as the pinch technique (PT) [t| to obtain a gauge independent WWV* Green function satisfying the 
requirements of gauge invariance and infrared finiteness. The PT exploits the fact that although the off-shell Green 
functions are gauge dependent, the S'-matrix elements to which they contribute are gauge independent. In this way, 
order by order in perturbation theory, one can construct a gauge independent Green function for an off-shell vertex 
by combining its contribution with all other Feynman diagrams that also contribute to a particular physical process, 
thereby getting rid of any £ dependent terms. In this respect, while the WWV* vertex does not represent a physical 
process by its own, it can contribute to a gauge invariant physical process such as e + e~ — > W + W~ scattering. By 
this method, a gauge invariant and gauge independent WWV* Green function can be constructed. The PT was 
also used in Ref. [T(| to calculate the one-loop contributions to the WWV* vertex from the unconstrained minimal 
supersymmetric standard model (MSSM). 

The PT is thus a diagrammatic method that allows one to remove any unphysical gauge dependent term from 
an off-shell n-point function at the level of the Feynman graphs contributing to a certain physical process, thereby 
yielding a well-behaved gauge-independent and gauge-invariant Green function. The PT may turn itself, however, 
into a somewhat cumbersome task. There is also the alternative of removing any gauge dependent term at the level 
of the generating functional. In this respect, although the background field method (BFM) [ll| renders a gauge 
invariant quantum action, no mechanism has been found yet that allows one to obtain both gauge-invariant and 
gauge-independent Green functions from the generating functional. Eventually, a connection between the PT and the 
BFM has been found it turns out that the gauge- invariant and gauge- independent Green function obtained via 
the PT is exactly reproduced if it is calculated through the BFM using the Feynman rules given in the Feynman-'t 
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Hooft gau ge ( £ = 1). Although such a correspondence, which so far remains a puzzle, was first established at the one- 
loop level |l3j , it has been shown that it persists at all orders of perturbation theory [l4| • Therefore, instead of using 
the PT, one can use the BFM Feynman-'t Hooft gauge (BFMFG) to calculate gauge-independent off-shell amplitudes. 
Along these lines, in Ref. [l[ we calculated the one-loop contributions to the WWV* vertex from the gauge sector 
of the so called minimal 331 model [HI, [l6|. In order to construct a well-behaved gauge independent Green function, 
rather than using the PT, we invoked an alternative method inspired in the BFM and BRST symmetry [T3|, which is 
well suited to study the sensitivity of the WWV* vertex to the virtual effects of the new gauge bosons predicted by 
331 models. We will follow a similar approach here to calculate the one-loop contributions from the 331 model with 
right-handed neutrinos [l8j . Although we will shortly describe the calculation scheme below, for more details we refer 
the reader to Ref. [lj. 

331 models have been the source of great interest lately. Recent studies within the framework of this class of 
theories have focused on neutrino physics [l9[ , Z' physics [2(| , Higgs boson physics [2l[ , bilepton physics P, @, [H, [H[ , 
supersymmetric extensions [24[ , dark matter [25| , CP violation [26[ , and theoretical aspects [13] • Models of this kind 
are based on the SUl(3) x Ux(1) gauge symmetry [l(| and are unique in the sense that they require that all 
the 3 fermion families be summed up in order to cancel anomalies, in contrast with other models in which anomaly 
cancelation is achieved family by family. As a consequence, 331 models require that the number of fermion families be 
a multiple of 3, the quark color number. This may suggest a solution to the family replication problem. Apart from 
this feature, 331 models are interesting as they predict a pair of massive gauge bosons arranged in a doublet of the 
electroweak group. While the minimal 331 model predicts a pair of singly-charged and a pair of doubly-charged 
y ±± gauge bosons, the model with right-handed neutrinos predicts a pair of neutral no self-conjugate gauge bosons 
Y° instead of the doubly-charged ones. These new gauge bosons are called bileptons since they carry two units of 
lepton number, thereby being responsible for lepton-number violating interactions [l8j . The neutral bilepton has 
been deemed a promising candidate in accelerator experiments since it may be the source of neutrino oscillations [28| . 
The reason why the effects of the bileptons on the WWV* vertex are worth studying is because their couplings to 
the SM gauge bosons are rather similar to the SM gauge boson couplings. Current bounds put the bilepton masses 
in the range of a few hundred GeVs [29| . which means that the bileptons may show up via their virtual effects 
in low-energy observables. This is an important reason to investigate the effect of these particles on the WWV* 
vertex. Furthermore, due to the spontaneous symmetry breaking (SSB) hierarchy, the splitting between the charged 
and neutral bilepton masses is bounded by m^, so the bilepton masses might be almost degenerate since they are 
expected to be heavier than mw- We will see below that this is particularly suited for our calculation method. 

The rest of the paper is organized as follows. In Sec. II we present a brief introduction to the 331 model with 
right-handed neutrinos. A survey of the quantization method is presented in Sec. Ill along with a detailed discussion 
on the gauge fixing procedure for the bileptons, whereas Sec. IV is devoted to the analytical results and the analysis. 
Finally, the conclusions are presented in Sec. V. 

II. THE 331 MODEL WITH RIGHT-HANDED NEUTRINOS 

The 331 model with right-handed neutrinos was first introduced in Ref. [181 ]. In a previouspaper, we worked in 
the context of this model and calculated the static electromagnetic properties of the W boson [7| . Also, the one-loop 
contributions to the static electromagnetic properties of the neutral no self-conjugate Y° boson were calculated in 
Ref. [23]. It was shown that the main contributions to the WW*y vertex arise from the gauge sector, i.e. from the 
bileptons, as the fermion sector does not contribute and the Higgs sector gives a negligibly contribution quite similar 
to that arisen in a THDM [7[. Therefore, although for completeness we will present a short description of the main 
features of the model, the following discussion will be mainly focused on the gauge sector as it is the one which is 
expected to give the dominant contributions to the WWV* vertex. Furthermore, we will see that our calculation 
scheme is suited to analyze the effects of the bileptons on the WWV* Green function. 

In the fermion sector of the 331 model with right-handed neutrinos, leptons of the same generation are arranged 
in left-handed triplets and right-handed singlets, the same being true for each quark generation. Apart from the 
introduction of right-handed neutrinos, there are three new quarks D\, D2, and T, which have the following electric 
charge Qd, = — l/3e and Qt — 2/3 e. In order to cancel the SUl(3) anomaly, two quark families transform as 
SUl(3) antitriplets and the remaining one as a triplet. At the first stage of SSB, when the SUl(3) x C//v(1) group is 
spontaneously broken into SUl(2) x LV(1), the new quarks get their masses and emerge as singlets of the electroweak 
group. Consequently, they cannot interact with the W boson at the tree level. It follows that there is no contribution 
from the new quarks to the WWV* vertex at the one-loop level. 

As far as the scalar sector is concerned, several Higgs sectors have been proposed in the literature to achieve the 
SSB in 331 models 13011. As for the 331 model with right-handed neutrinos, it requires the simplest Higgs sector of 
this class of theories [lg, i.e. only three triplets of SUl(3) are required to reproduce the SM physics at the Fermi 
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scale: 



X=(*o)~ (1,3,-1/3), P= (1-3,2/3), ^7= (^o) ~ (1,3,-1/3), (1) 

where $ y = (Gy , Gy ) , $J = 0~,p *) and $2, = (ry°*,?7 + ) are SE/ L (2) x J7y(l) doublets with hypercharge -1, 1 
and -1 respectively. The vacuum expectation values (VEV) are < x > T = (0, 0, w/V2), < p > T = (0,u/V2,0), and 
< rj > T = (v/y/2, 0, 0). The triplet x breaks the SUl(3) x Un(X) group into SUl(2) x t/y(l) at the u> scale, whereas 
p and r\ are meant to break SUl(2) x t/y(l) into U e (l). 

In the gauge sector, the model predicts the existence of five new gauge bosons: two singly charged Y^, two neutral 
no self-conjugate Y°, and a neutral self-conjugate Z'. All these new gauge bosons and the new quarks as well acquire 
their masses at the w scale. At this stage of SSB, the Z' boson emerges as a singlet of the electroweak group and so no 
interaction between the Z' boson and the W boson arises. However, in the following stage of SSB, at the Fermi scale, 
the Z' boson couples with the W boson via the Z — Z 1 mixing angle. As a consequence, the Z' contribution to the 
WWV* vertex will be proportional to the square of the Z — Z 1 mixing angle, which is highly suppressed according to 
the most recent experimental bounds [U, [32j- We will thus ignore this contribution in this work. In sharp contrast, 
the bileptons arise as a doublet of the electroweak group at the w scale, which means that they couple to the SM 
gauge bosons in a rather peculiar way: due to the fact that the SUl(2) group is completely embedded in SUl(3), the 
bileptons couplings to the SM gauge bosons have similar strength and Lorentz structure than those of the SM gauge 
boson couplings themselves. These new interactions, which arise solely from the Yang-Mills sector, do not involve any 
mixing angle and are entirely determined by the SUl{2) coupling constant and the weak angle Qw 

In the first stage of SSB, the bilepton masses are degenerate: myo = m Y ± = my = gw/2. However, when 
SUl(2) x Uy(X) is broken down to U e (l), the bilepton masses degeneration is also broken. Once the Higgs kinetic- 
energy sector is diagonalized, there emerge the following mass-eigenstate fields: 
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with masses m Y0 — g 2 {w 2 + u 2 )/4, m Y± — g 2 (w 2 + v 2 )/A, and = g 2 {u 2 + u 2 )/4. The remaining three gauge 
fields A^, A^,and define the self-conjugate mass eigenstates. From the above expressions, it is straightforward to 
obtain the following upper bound on the bilepton masses splitting [l8l |: 

\m Y o - Wy± I < mly. (5) 

Therefore, the bilepton masses become nearly degenerate when they are much larger than mw- While in the minimal 
331 model an upper bound on the bilepton masses of the order of 1 TeV can be derived from the theoretical constraint 
sin 2 9\\r < 1/4 [l6t I321. [33| . which is obtained from matching the gauge couplings constants at the SUl(3) x Ux(1) 
breaking scale, in the version with right-handed neutrinos the corresponding bound is highly relaxed as the theoretical 
constraint is sin 2 0w < 3/4 

The Yang-Mills sector induces all the couplings we need to compute the bilepton contribution to the WWV* vertex, 
so we will devote particular attention to it. The respective Lagrangian can be written as 



C. Y m = -\f; v F^-\n iiv N^, (6) 

where F* — d^A^ — d w A a ^ + f^A^Af, and — d^N v — d w N^ : being f abc the structure constants associated with 
SUl(3)- When the mass eigenstate fields are introduced, the above Lagrangian can be split into three SUl(2) x Uy (1)- 
invariant terms [l8| : 

r _ r SM , r SM-NP , r NP ( ~\ 
J-YM — l-YM "+" %M "+" L -YM- \ ' ) 

While C y A m stands for the usual SM Yang-Mills Lagrangian, C Y A ^f NP comprises the interactions between the new 
gauge bosons and the SM ones. The latter Lagrangian is the only one relevant for the present discussion and we will 
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get back to it later. As for the last term, CyMi ^ induces the couplings of the Z' boson to the bileptons. At the 
one-loop level there are no contributions to the WWV* vertex induced by this Lagrangian. 

In order to calculate the one-loop correction to the WWV* vertex, we will introduce a SUl(2) x Uy (l)-covariant 
gauge-fixing procedure for the bileptons. This will be discussed in the following section. 

III. A SU L (2) x U Y (l) COVARIANT GAUGE-FIXING PROCEDURE FOR THE BILEPTONS 

A. Overview of the quantization method 

We now turn to present an overview of the quantization method used to obtain a gauge invariant and gauge 
independent Green function for the WWV* vertex in the 331 model with right-handed neutrinos. This method, 
which is inspired in the BFM and BRST formalism, was comprehensively discussed in our previous work [l[ and we 
refrain from presenting a detailed discussion here. 

The BFM is a powerful tool that allows one to construct a gauge invariant quantum action out of which 
gauge invariant Green functions can be obtained that are free of pathologies and satisfy simple Ward identities. In 
a conventional quantization formalism, all the fields appearing in the action are quantized. In the BFM, the gauge 
fields are decomposed into a quantum field and a classical (background) field. While the quantum fields are integrated 
out, the background fields are treated as sources. As a result, the quantum fields can only appear as internal lines in 
loop diagrams, whereas the background fields appear as external lines. In principle, it is necessary to gauge-fix both 
the quantum and the classical fields in order to define S'-matrix elements, but it is enough to gauge-fix the former to 
obtain off-shell Green functions. In this respect, while gauge invariance with respect to the quantum fields is broken 
when they are gauge-fixed, this process leaves unaltered the gauge invariance with respect to the classical fields. The 
resulting Green functions will be gauge invariant but they will still maintain the dependence on the gauge parameter 
that characterizes the gauge-fixing scheme used for the quantum fields. However, we can exploit the connection 
between the BFM and the PT to obtain a gauge invariant and gauge independent Green function. We just need to 
use the BFM Feynman rules given in the Feynman-'t Hooft gauge. We will see below that our quantization method 
incorporates some features of the BFM. 

Although gauge invariance with respect to the full theory is broken when a conventional quantization scheme is 
applied, one can still preserve the gauge invariance under a subgroup of the theory. This approach is well suited when 
the virtual effects of the gauge fields associated with the subgroup are expected to be considerably small. In this 
context, following the philosophy of the effective Lagrangian approach, where a SUl(2) x LV(1) effective Lagrangian is 
constructed out of the light (SM) fields to assess the effects of the heavy fields on a low-energy observable, it would be 
convenient to analyze the virtual effects of the bileptons on the WWV* Green function in a SUl(2) x LV(l)-covariant 
manner. A quantization scheme for the bileptons would only be required since the SM gauge fields would only 
appear as external lines. Thus, a SUl{2) x Uy (l)-invariant effective Lagrangian can be constructed by introducing a 
SUl{2) x t/y (l)-covariant gauge-fixing procedure for the bileptons, which then can be integrated out. The gauge-fixing 
procedure must involve the SUl{2) x Uy (l)-covariant derivative given in the representation in which the heavy fields 
transform under this group, which is the reason why such gauges are known as nonlinear or covariant. This class of 
gauges have proved to be very useful for the calculation of radiative corrections in the SM and beyond [lij . 

Since we want to introduce a nonlinear gauge fixing procedure for the bileptons, we need to be careful as the 
difficulties to implement the Faddeev-Popov method (FPM) 35] in such a case have long been known. More specifically, 
it is known that renormalizability becomes ruined when the FPM is applied to a nonlinear gauge. We can invoke 
instead the BRST formalism [17| to construct the most general renormalizable nonlinear gauge-fixing term [36j |. We 
will thus introduce a SUl(2) x Uy (l)-covariant gauge-fixing term for the bileptons, which will lead to an invariant 
quantum action out of which a gauge invariant WWV* Green function will be obtained. Invoking the connection 
between the PT and the BFM, the use of the Feynman-'t Hooft Feynman rules will render a gauge independent Green 
function. 

Finally, we would like to point out that both our quantization method and the BFM are meant to construct gauge- 
invariant quantum actions. The main difference resides in the fact that while the BFM allows one to analyze any 
new physics effects regardless the energy scale, ours is only appropriate to study heavy physics effects on low-energy 
(SM) Green functions. This stems from the fact that while in the BFM gauge-invariance would be preserved with 
respect to the gauge group of a complete theory, in our quantization method there is only invariance under a subgroup 
of such a theory. In the case of the present paper, our quantization method preserves gauge invariance under the 
SUl(2) x C/y(l) group rather than SUl(3) X Ujy(l). While the calculation will be greatly simplified because of 
electroweak invariance, the price to be paid is that our result will only be approximate as it will only account for the 
bilepton effects on the WWV* vertex at the w scale, when the bilepton masses are still degenerate. We will see below 
that this is indeed a good approximation. 
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B. Gauge-fixed Lagrangian and Feynman rules 

The procedure to gauge- fix the Yang- Mills Lagrangian of the 331 model with right-handed neutrinos is similar to 
that described in the case of the minimal 331 model. So, we refrain from presenting a detailed discussion here and 
refer the reader to Ref. [l[. 

After introducing the most general action for a Yang-Mills system consistent with BRST symmetry and renor- 
malization theory [37| . we integrate out the auxiliary fields to obtain an action defined by the following gauge-fixed 
SUl(2) x Uy (l)-invariant Lagrangian 

CbRST = C-YM + C-GF + £-FP, (8) 

where Cgf and Cfp are the gauge-fixing term and the ghost sector Lagrangian, respectively. As for the Yang-Mills 
Lagrangian Cym, which is given in Eq. 0, for our calculation we only need the £y-^ _Arp term as it is the only one 
that induces the interactions between the bileptons and the SM gauge bosons. It can be expressed as 



J0 Y M~ NP = -\ (D»Y» - D v Yrf {D»Y V - D V Y») - iY} (gF^ + g'B^) Y v 

_ 4s 2 w / f yu _ D „ _ yv _ DU t \ 

2 cw ^ ' 

where YJ = (Y®*,Y+) is a doublet of the electroweak group with hypercharge — 1, and Z? M = — igA^ + ig'Tl^ 
is the electroweak covariant derivative. We have also introduced the definitions F^ v = a l Ff ll/ /2, A M = a t A] i /2, and 
Bp = YB^/2, with <j % the Pauli matrices. 

As for the gauge- fixing term Cgf and the ghost sector Lagrangian Cfp, they can be written as 

£ G F = -^rr, (io) 

and 

C F p = -C a {5f a ) - jf sbo f s C b C e - ^f abe f cde C s C b C i C e , (11) 

where f a are the gauge-fixing functions and C a stands for the antighost fields. In addition, C a are the ghost fields 
associated with the A" fields, f abc are the SUl(3) structure constants, and £ is the gauge parameter. 

According to Ref. [l|, the most general SUl(2) x Uy (l)-covariant gauge-fixing functions f a consistent with renor- 
malization theory is given by: 

f = (S u %-gf u Ai)A^ b -^f bs x' i X b X, a = 4,5,6,7; i= 1,2,3,8. (12) 

We can now insert this expression into the gauge-fixed Lagrangian. Apart from analyzing the dynamics induced by 
each term of the gauge- fixed Lagrangian, we would like to put special emphasis on the covariance under SU i (2) x Uy (1) • 
The covariant structure of the gauge-fixing term becomes manifest when the mass eigenstates f Y = (/ 4 — if 5 ) 

and fy = 4= (/ 6 — if 7 ) are introduced in an SUl(2) x Uy(1) doublet 



We will now decompose the gauge-fixing Lagrangian into three terms in order to analyze the dynamics it induces: 

Cgf = Cgfi + Cqf-2 + ■ ■ ■ , (14) 
where each term is separately SUl{2) x Uy(1) invariant: 

Cgfi = -C\D^)\D V Y V ) - ^( X ' * X '°)^ Y <S> Y ), (15) 



C G F2 = ^ (x°*(D^)^y x *UWj) , (16) 
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whereas the third term, denoted by . . . , involves the Z' boson and is irrelevant here. 

We note that the first term of Cqfi defines the bilepton propagators and gives new contributions to the trilinear 
and quartic couplings involving the bileptons and the SM gauge bosons, which originally arise from the CyM~ N 
Lagrangian. After some algebra, we are lead to the Feynman rules for these modified couplings, which are shown in 
Fig. [TJ The Lorentz structure associated with the trilinear couplings is given by 

r QM „(fc, ki,k 2 ) = (k 2 - h) a g fiv + (k-k 2 - £,~ l k\) g av - (k - k x - £~ 1 k 2 ) u g a n, (17) 

whereas those of the quartic couplings are 

r SftL yr = 9apg^ - Zg^gpn + (l + C 1 ) g a ^gp Vl (18) 

and 

= (Qy- + QY°)g a p9v, + [(1 + r 1 ) Qy- - 2Qy~-] g afl g^ + [(1 + r 1 ) Qy- - 2Qy«] g av g^. (19) 

The reason why the trilinear couplings WYY and VYY have the same Lorentz structure is a consequence of 
SUl(2) x UyiX) invariance. This is also reflected in the fact that a simple Ward identity is satisfied by these vertices: 

k a Y aliV {kM,k2) = nll Y \k 2 )-n%J(k 1 ), (20) 

where H YY (k) is the two-point vertex function 

n£» = (-* + m Y ) 9tw - (r 1 - 1) k,K. (21) 





FIG. 1: Feynman rules for the trilinear and quartic vertices involving the bileptons and SM gauge fields in the SUl{2) x 
covariant _R^-gauge. ev = e, Qy- = — 1, and Qyo = for V = 7, whereas ev = g/(2cw), Qy- = 2s^- — 1, and Qy = 1 for 
V = Z. In addition, e w — g/y/2. See Eqs. (ri7|) - (fl9jl for the Lorentz structures. 

Note that the unphysical masses for the pseudo-Goldstone bosons associated with the bileptons are determined by 
the scalar part of £gfi, which also modifies various unphysical couplings arising from the Higgs potential. As for 
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£-GF2, it allows one to remove various unphysical vertices that arise from the \ kinetic-energy sector Ck x , which can 
be written as: 



where each term is SUl(2) x Uy (l)-invariant by itself: 



(22) 



(23) 



C 



K X 2 = 



X°*YUD^ y ) + <P Y Y^x a - H 



(24) 



Once again, the terms denoted by . . . are irrelevant for our calculation as they involve the Z' boson. 

From these expressions it is clear that Ck x i induces the couplings of the pseudo-Goldstone bosons to the SM gauge 
bosons. The respective Feynman rules for these coupling can be extracted straightforwardly and are presented in Fig. 
[3J As for Ck X 2, it induces the bilinear terms Y^ *Gy* '° and Y^G Y , together with the unphysical trilinear couplings 
Y®> *W~ ,+ G Y ' ~, Y+'~W~' + G Y a * , and some quartic vertices irrelevant for our calculation. It turns out that all these 
couplings exactly cancel when Lk X 2 and Cqf2 are combined: 



C KX 2 + £ G F2 = ie w [x °*^(y^$ y ) + tfyY^x - H.c. 



(25) 



where . . . stands for irrelevant surface terms. This is the reason why this gauge-fixing procedure will simplify con- 
siderably our calculation as we can get rid of several Feynman diagrams involving the couplings Y®*°*W~' + G Y '~ and 



Y+'-W-'+G /*. 



Finally, we would like to show the covariant structure of the ghost sector and extract the Feynman rules necessary 
for our calculation. We introduce the mass eigenstates in SUl(2) x E/y(l) doublets: 



Cy 

cz 



Cy 



Cy 

c; 



(26) 



where the mass eigenstates are defined as C Y — -^(C 4 — iC 5 ) and C Y — "^(C 6 — iC r ), and similar expressions for 
the antighost field Cy. The Cfp Lagrangian can thus be written as follows: 

C FP = (D^C Y )HD^C Y ) + 9 - [{YyY»){C Y a l C Y ) + 3(Y^)(C Y C Y ) - A{Y^C Y ){Y^C Y ) 



where 



+ ^Y}M C D»C Y + ^YlMcCY f 



V2 
i\/2\ 



X°*X°C y Cy + x°$y m cCy - (C y $y)($ y CV) 



(M C C Y + M c C y )Hd^) - (D^){M C C Y + M C C Y ) 



$ y (M c C y + M c Cy)x° + X°*(M C C Y + M c C Y y$Y 



H.c. 



j?abc jcde ^ja ^jb ^jd ^je 



(27) 



M c 



^(c 3 +v3c*: 
kc^c 2 ) 



V2 



^-ic*) 

-L(C 3 -V3C*) 



(28) 



Mt 



(29) 



where i = 1,2,3,8, T) l -> = S^d^ — gf l ^ a A ( t is the covariant derivative given in the adjoint representation of SUl(3), 
and Mc is obtained from Mc after replacing the ghost fields by the antighost fields. Under a SUl(2) x U y (1) unitary 
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y y 




FIG. 2: Feynman rules for the trilinear and quartic vertices involving SM gauge fields and scalar unphysical particles (pseudo- 
Goldstone bosons and ghosts) in the SUl(2) X Uy (1 )-covariant i?j-gauge. In this gauge, the W and V couplings to pseudo- 
Goldstone bosons and ghosts coincide. 



transformation U, Mq transforms as Mq — > UMqU^ . The same is true for Mq and M-c- It is thus clear that Cpp 
is SUl(2) x £/y(l) invariant. 

From Cfp, it is straightforward to show that the Feynman rules for the trilinear and quartic couplings involving 
the ghost fields and the SM gauge bosons have the same Lorentz structure than those involving the pseudo-Goldstone 
bosons and the SM gauge bosons, which stems from the fact that each sector is SUl(2) x Uy(1) invariant by itself. 
The Feynman rules for the ghost (antighost) fields are summarized in Fig. [2] together with the Feynman rules for the 
pseudo-Goldstone bosons. We see that the trilinear vertices WS>S and VS'S, with S standing for a commutative 
(pseudo-Goldstone boson) or anticommutative (ghost) charged scalar, satisfy simple Ward identities: 

k a Tl sts =IL s ' s '(k 2 )-U ss (k 1 ), (30) 

where r^ 5 * 5, = (fti — k 2 ) a and IL ss (ki) stands for the two-point vertex functions n(fc) = fc 2 — ^m Y . 



IV. ANALYTICAL RESULTS AND DISCUSSION 

A. A gauge invariant and gauge independent Green function for the WWV* vertex 

The most general transverse CP-even vertex function for the W+{p — q)Wg(—p — q)V*(q) coupling has the form 
0,11 

r I/3 M = - i 9v{A[2p Pi g a/} + 4(qgff aM - g a g/3 M )] + 2An v (q fJ g atl - q a g 0fl ) + 4A ^ y (p^q a q fj - \q 2 p^g a p)\- (31) 

The longitudinal terms were ignored as they become negligibly small when V* couples to a light fermion current, 
as in e~e + — > WW scattering. The SM tree-level values are A = 1 and An = AQ = 0. While An, and AQ are 
ultraviolet finite at the one- loop level, A is divergent and requires renormalization. The coefficients Ak 7 and AQ 7 , 
obtained from the WW7 vertex function, determine the static electromagnetic properties of the W boson, namely, 
its magnetic dipole moment fiw and its electric quadrupole moment Qw- 

From the Feynman rules presented in Figs. 1 and 2, we can construct all the Feynman diagrams contributing to the 
WWV vertex at the one-loop level, which we show in Fig. 3, and extract the Any and AQv coefficients. Although 
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there is a similar set of Feynman diagrams for the pseudo-Goldstone bosons and the ghost fields, the unphysical 
particles only contribute to Eq. (f3"Tj) via the triangle diagram. Moreover, it turns out that the ghost (antighost) 
contribution is minus twice that of the pseudo-Goldstone bosons, which is due to the separate SUl(2) x Uy(1) 
invariance of these sectors. We have evaluated the loop amplitudes in the Feynman- 't Hooft gauge via the Passarino- 
Veltman reduction method [3^|. We have verified that the bilepton contribution can be written exactly as Eq. (|3ip . 



which means the WWV* Green function satisfies the Ward identity T^^^g^ = 0. It turns out that the form factors 



are the same regardless of V, which means that, just as occurs at the tree-level, and rfg only differ by 



^d rf, 

the gv coefficient. Thus SUl(2) x Uy(1) invariance is preserved at one-loop level. The gv coefficient is given by 
gv = ev(Qyo — Qy- )■> which, after inserting the charge values, gives gv = e for V = 7 and gv = cwg for V = Z. 

We now introduce the following short-hand notation Q = 2q, Q 2 = Q 2 /m 2 v , iy = m Y /m 2 v , (3(Q 2 ) = 6a/(4 — Q 2 ) 3 , 
a = g 2 /967r 2 , and express the A^y and AQv coefficients in terms of two- and three-point Passarino-Veltman scalar 
functions: 



AQ V = 2/3(Q 2 ) 



(6x Y G(Q 2 ) + 1) Q 4 - 2 (2 (2x Y + 1) F{m 2 w ,Q 2 ) + 2x Y F(0, Q 2 ) + 6 (3x Y - 1) G(Q 2 ) + 3) Q 2 



+ 4 ((8a* - 5) F(m 2 w ,Q 2 ) + 8x Y F(0, Q 2 ) + 6 (3x Y - 1) G(Q 2 ) + 5) 



8 (3F(m 2 w , Q 2 ) - 8x Y F(0, Q 2 ) + 3 (1 - Ax Y ) G(Q 2 ) - 6) 

Q 



(32) 



Ak v = (3{Q 2 )y{F{m 2 Wl Q 2 ) - Ax Y F{m 2 w ,Q) - Q (l + 3x Y ) G(Q 2 ) - 3) Q 4 + 2 (2 + 16x Y F(m 2 w , 0) 

+ (13 + I6^)F(m^,g 2 )-3(l-3xy)G(Q 2 ))Q 2 + 32(l-2 (F(m 2 w ,0) + 2 F(m 2 w ,Q 2 ))) x Y ], (33) 

where F(m 2 ,Q 2 ) = B (m 2 ,m Y ,m Y ) - B Q (Q 2 , m\ , m\) and G(Q 2 ) = C (Q 2 ,m 2 Vl m 2 v ,m\,m\,m\), with the 
scalar two-point Bq and three-point Co functions given in the usual notation. From here, it is clear that ultraviolet 
divergences cancel up in the W form factors. 

Before the numerical evaluation, it is interesting to analyze the static electromagnetic form factors, i.e. the scenario 
when Q 2 = 0. Note that Q 2 and F(0, Q 2 ) vanish when Q 2 = 0, whereas G(0) = (2 - F{m 2 w , 0))/(l - ^x Y ) [40j. After 
taking the Q 2 — * limit in the last two terms of Eq. ([32|) . AQ 7 can be written as 



AQ n 



§m 2 w (\ - 4x Y ) 



(-21 + A8x Y + 3 (I + (10 - 32x Y ) x Y F{m 2 Wl 0))) 



+ (82* -3) 



dB Q (Q 2 ,m Y ,m Y ) 



dQ 2 



+ imly (1 - Ax Y 



dC {Q 2 , m^y, mly, m Y , m Y , m Y ) 



Q 2 =0 



dQ 2 



, (34) 



Q 2 =0, 



where L'Hopital rule has been used for the indeterminate limit. Since any three- and two-point scalar function and 
their derivatives can be written as a combination of two-point functions [40( , we can write 



dB Q (Q 2 1 m Y , m 



1 



dQ 2 



3 2 =0 



6m Y ' 



and 



dCQ{Q 2 ,m 2 v ,m 2 v ,m\,rri Y ,m Y ) , 



dQ 2 



? 2 =o 



(5 — 6x Y ) iyF(m^, 0) — 5x Y 
(1 - Ax Y ) 2 ' 



We thus obtain 



AQ 7 = 



4a 



1 - 4x Y 



(-1 + (1 - 3 (2x Y - 1) (B {m 2 w , m 2 Yl to 2 ,) - B (0, m 2 Yl to 2 ,))) x y ) 



(35) 



(36) 



(37) 



whereas Ak 7 can be obtained immediately from (|33 



Ak^ 



3« 



(l - 6x Y (Boimly ,m Y ,m Y ) - B o (0,m Y ,m Y ))) . 



(38) 
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Y~ Y~ 



FIG. 3: Feynman diagrams for the WWV vertex in the SUlC2) x Uy (l)-covariant gauge. The pseudo-Goldstone bosons and 
ghosts contribute through an identical set of diagrams, but only the triangle ones give a nonvanishing contribution to Any and 
AQv 

In the limit of very large xy (my 3> m\y) we obtain 

B (™^,TOy,my) - B Q (0, m\,m\) ~ — - , (39) 

OXy 

which can be used to show that both AQ 7 and Ak 7 vanish in the large bilepton mass limit, i.e. these coefficients are 
of decoupling nature. Furthermore, it is evident that AQ 7 decreases more rapidly than Ak 7 as my increases. 

We emphasize that the above expressions account for the bilepton effects on the WWV* vertex at the w scale, 
when the bilepton masses are degenerate and equal to my . These results are thus approximate but can give a realistic 
estimate of this class of effects. We now proceed to the numerical evaluation and analysis. 

B. Numerical evaluation and discussion 

We now would like to analyze the behavior of AQv and Any as functions of Q 2 and the bilepton mass, which are 
the only free parameters which they depend on. 

As far as Q 2 is concerned, we will focus on the values that can be at the reach of the future planned linear colliders 
as this vertex has the chance of being studied through e~e + — > W~W + scattering. We thus consider the range 200 
GeV < y/Q 2 < 2 TeV 0. 

As for the bilepton mass, it is convenient to give a short account on the existing bounds from both the theoretical 
and experimental sides. First of all, we already mentioned that in the 331 model with right-handed neutrinos the 
matching of the gauge couplings constants at the SUl{3) X Ux(l) breaking scale requires that sin 2 6y/ < 3/4 18], 
from which an upper bound on the bilepton mass can be derived after considering radiative corrections to sin 9w 
and the most recent experimental data [42j]. The respective bounds is somewhat weak and contrasts with the case 
of the minimal 331 model, which requires that my < 1 TeV 0, H2, HH f° r consistency with the theoretical bound 
sin 2 9\v < 1/4. Therefore, whereas the bileptons predicted by the minimal 331 model could be searched for via 
collider experiments in a near future, thereby confirming or ruling out the model, the 331 model with right-handed 
neutrinos would still leave open the door. We have also mentioned that because of the symmetry breaking hierarchy, 
the bilepton mass splitting is bounded by | rriy^ — m^ 1 < ttiw- This means that myo and my± cannot be arbitrarily 
different. In fact, the heavier the bilepton masses, the closer its degeneracy. Our approximation of mass degeneracy 
has thus much sense. We now would like to examine the lower bounds on the bilepton masses. In Ref. [2!| it 
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was argued that the data from neutrino neutral current elastic scattering give a lower bound on the mass of the 
new neutral gauge boson mz> in the range of 300 GeV, which along with the symmetry-breaking hierarchy yield 
royi ~ m,Y ~ 0.72 mz' > 220 GeV. A similar bound was obtained in Ref. [28| from the observed limit on the wrong 



muon decay R = r(/j, 
These lower bounds on 



) < 1.2%, which leads to m Y ± > 230 ± 17 GeV at 90% C.L. 



are in agreement with that obtained from the latest BNL measurement on the muon 



anomaly [29j. It is then reasonable to consider the range 200 GeV < my a < 2000 GeV for our numerical analysis. 
This will allow us to assess the behavior of the form factors and get an estimate of their order of magnitude. A word 
of caution is in order here: it has been pointed out that some of the existin g bo unds on the bilepton masses are too 
model dependent and thus can be relaxed by considering extra assumptions [43J , so a relatively light bilepton cannot 
be ruled out yet. 

In Fig. 5] we show the W form factors Any and AQv vs the center-of-mass energy Ecm = \J Q 2 , i.e. the momentum 
carried by the virtual V boson, for the bilepton mass values my — 200, 600, and 1000 GeV. On the other hand, the 
dependence of the W form factors on the bilepton mass is shown in Fig. where we have plotted Any and AQv vs 
the bilepton mass for Ecm = 500 and 1000 GeV. 



2300 3000 




FIG. 4: The W form factors AQv (upper plots) and Akv (lower plots) vs the center-of-mass energy Ecm = a/ Q 2 for various 
values of the bilepton mass. The form factors are in units of 10~ 6 and Ecm is in units of GeV. The solid (dashed) lines 
represent the real (imaginary) part of the form factors. 



As far as the energy dependence of the form factors, from Fig. 0]it is clear that both AQy and Any are of the order 
of about 10~ 4 for a relatively light bilepton with a mass of 200 GeV and a center-of-mass energy around the energy 
threshold Ecm = 2my, where the form factors reach their extremum values and develop an imaginary part. Below 
this threshold, the form factors are purely real, which reflects the fact that the bileptons that couple to the V boson 
are necessarily virtual, whereas the direct production of a bilepton pair becomes possible for energies Ecm > 2my. 
Although the direct production of a bilepton pair would be preferred over the search for their virtual effects, the latter 
would be useful for a high precision test of the WWV vertex. We can also see that the form factors decrease by one 
order of magnitude when my = 500 GeV and by two orders of magnitude when my = 1000 GeV, which is depicted 
in Fig. [5j For very high energies, the form factors have a good behavior as they approach zero asymptotically after 
reaching an extremum value above Ecm = 2my. Therefore, unitarity is respected, which stems from the fact that 
the WWV Green function that we have obtained is gauge invariant. In Fig. Q] we can also observe that, except for 
the reversed sign, the curves for AQv and Akv are very similar. For a relatively light bilepton both form factors are 
of about the same order of magnitude, but the magnitude of AQv decreases more quickly as m Y increases. 

It is interesting to compare our results with those obtained in the SM and some of its extensions. As far as the 
SM is concerned, the gauge boson contribution to the AQv form factor is of the order of 10~ 3 — 10~ 4 for Ecm in 
the range 200 — 1000 GeV, whereas AQv is about one order of magnitude below 0, H[. These contributions are of 
the same order of magnitude than those of the fermion and scalar sectors of the theory 3[ . As far as supersymmetric 
models are concerned, the total contribution from the unconstrained MSSM is about the same order of magnitude or 
larger than the SM total contribution for some set of the values of the parameters of the model [l(| . Our results in 
the 331 model with right-handed neutrinos and relatively light bileptons, with a mass around 200 GeV, is of the same 
order of magnitude than the SM contribution. These results are similar to those obtained in the minimal 331 model 



12 




[l| and other weakly coupled theories. Unless a very high precision is achieved in future particle colliders, it would be 
very hard to discriminate the contributions form different models to the WWV vertex. However, testing this vertex 
would still be useful for probing some particular model once it has been confirmed by experimental data. 



V. FINAL REMARKS 



We have used a nonconventional quantization method inspired on the BFM and BRST symmetry to analyze 
the effects of the new gauge bosons (bileptons) predicted by the 331 model with right-handed neutrinos on the 
of-shell WWV* vertex. Hopefully, this class of effects might be searched for at a future linear collider through 
e~e + — * V — > W~W + scattering. In particular, it has been pointed out that CLIC would reach a sensitivity of 
about 10~ 4 in the measurement of the so called form factors, AQy and Any, characterizing the WWV* vertex [4l| . 
Following the philosophy of the BFM, our method considers the bileptons as quantum fields and the SM gauge bosons 
as classical fields. A nonlinear SUl(2) x UyiX) invariant gauge- fixing is then introduced for the bileptons, which are 
integrated out in the generating functional. The result is a quantum action defined by a SUi,{2) x UyiX) invariant 
Lagrangian 1 out of which a gauge invariant and gauge independent Green function can be obtained. To this end, 
we made use of the link between the diagrammatic method known as the PT and the BFM, which establishes that 
the gauge invariant and gauge independent Green function obtained via the PT coincides with that obtained via the 
BFM as long as the Feynman-'t Hooft gauge Feynman rules are used when calculating the latter. We emphasize 
that our method is only approximate as the quantum action is only invariant under SUl(2) X Uy(l) rather than 
SUl(3) x Un(1). This method is suited to analyze the bilepton effects at the SUl(3) x Un(1) breaking scale, when 



1 Actually, this scheme is quite analogous to the effective Lagrangian approach, in which an SUl(2) x £/y(l) effective Lagrangian is 
constructed out of the SM fields to analyze the virtual effects of the heavy fields, which have been integrated out. 
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their masses are still degenerate. The advantages of our calculation scheme are twofold: the introduced nonlinear 
gauge-fixing term for the bileptons allows one to remove several unphysical vertices, which in turn allows one to get 
rid of several Feynman diagrams; on the other hand, preserving the electroweak invariance turn the calculation into a 
simple task as each sector of the theory gives a gauge invariant contribution satisfying simple Ward identities by its 
own. Once our method was applied, we obtained the form factors AQv and Any, which were analyzed for several 
values of the bilepton mass and convenient values of the momentum carried by the virtual V boson. It was found 
that the bilepton effects on the WWV* vertex are of the same order of magnitude than the SM and the minimal 
331 model contributions, provided that the bilepton mass is of the order of a few hundred of GeVs. For very heavy 
bileptons, the respective contribution to the WWV* is negligibly small. This indicates that it will be very hard to 
discriminate between different class of effects on the WWV* vertex arising from distinct models. However, if a high 
precision is achieved in future linear colliders, the WWV* vertex might serve as a probe to some particular model by 
then confirmed by other experimental data. The good behavior of the form factors at high energies was also discussed, 
for it is a reflect of the gauge invariance and gauge independence of the WWV* Green function obtained via our 
quantization method. 
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